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Abstract 

We derive the exact rate of pair production of oppositely charged scalar particles by 
a smooth potential step </>(x) oc tanh/cz in three dimensions. As a check we recover from 
this the known results for an infinitely sharp step as well as for a uniform electric field. 



1 Introduction 

Spontaneous pair production of oppositely charged particles was first discussed by Heisenberg 
CN| ■ and Euler [1] by making use of the Dirac picture of the vacuum. It also provided the solution 
of the Klein paradox of relativistic quantum mechanics [2] for the scattering on a potential 
step. Within second quantized field theory it was first discussed by Hund [3], as a precursor 
of the later famous calculation of Schwinger [4], whose generalization to gravitational fields 
was made by Hawking [5]. Schwinger calculated the one- loop effective action of QED in a 
constant electromagnetic field by the proper time method. The imaginary part of this effective 
0> action yields directly the probability of pair production from the vacuum. The result confirms 
the fact anticipated by Sauter [6] twenty years earlier that particles can pass through strong 
repulsive potentials without the exponential damping expected in quantum tunneling processes. 
^ ■ A detailed review and the relevant references can be found in Refs. [7] and [8]. 

Another method to describe the problem of pair production is due to [9,10]. It is based on the 
use of causal Green functions [11] to relate the particle production directly to the scattering 
process. It can be shown that the rate of pair production can be expressed as an ordinary 
energy-momentum integral over the logarithm of the reflection coefficient. This formula allows 
one to connect in a most transparent and efficient way the one-particle Dirac approach in 
which the Klein paradox first appeared with the second quantized field theory in which the 
problem was first solved satisfactorily. It has been used for many other developments including 
semiclassical approximations [12-14]. In the same papers [9,10], also the scattering and pair 
production processes by the Sauter potential 0(x) oc tanh kz were analyzed in detail. The exact 
solutions of the both Dirac and Klein-Gordon equations were found in this potential, and the 
causal Green functions were constructed to define the pair production rates of fermions and 
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bosons. However, the calculations were never carried on to derive an exact expression for the 
pair production rate. It is the purpose of this paper to complete this gap. 



2 Barrier scattering 

Consider a relativistic scalar particle of charge e > and mass m moving in an external 
electromagnetic potential A^(z) = (<f)(z),0) corresponding to nonuniform electric field along 
z-direction with the strength E{z) = —d(f>(z)/dz. Its kinetic energy and momentum are 



Po(z) = p -e</>(z), p 3 (z) = yJP$(z) -(pl + m 2 ), p\ = p\ + p\. (1) 

In the transverse direction the particle propagates freely as a plane wave exp[i(p_|_-x^ — p Q x )/h]. 
Thus we represent the solution of Klein-Gordon equation in the potential 0( product 
of this plane wave with the z-dependent wave function ip(z) satisfying the Schrodinger-like 
equation (c = 1): 

r(z) + jp [(po - e<P(z)) 2 - (p\ + m 2 )] ftz) = . (2) 

We shall consider the smooth step potential of the Sauter type 

e<j)(z) = v tanh kz , (3) 

with v, k > 0, where the sharp step potential is recovered in the limit k — > oo , while the limit 
k — > with fixed vk reproduces the linear potential due to a constant electric field. 

In the potential (3), Eq. (2) can be solvable exactly [10]. The solution describes the barrier 
scattering of a particle impinging the left with asymptotic boundary conditions 



1>(z) 



A l e*Pa )z ' h + A 2 e~ ip B )z ' h , z -> -oo 



B e ip z~ )z l % , z -»• +oo 



(4) 



where Ai, A 2 and B are normalization constants. Let us introduce the initial and final values 
of the particle energy far to the left and to the right of the potential (3): 

P!? ) = Po(z)\ z ^ 00 =po±v, (5) 
and the corresponding momenta 

= p 3 WU TOO = vW T)2 -(pi + ™ 2 )- (6) 

In the scattering process, the momenta p^ in Eq. (6) are real thus restricting the asymptotic 
energies to Pq T ^ > \jp\ + m 2 . 

To solve Eq. (2) with the asymptotic conditions (4) we set 

p^ = lUkix , p { 3 +) = 2hkv , (7) 
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where 



2 2 v Po 
a — v — — . 

^ h 2 k 2 

We further replace z by the dimensionless variable 

C = - exp(-2kz) , 

running from — oo to 0. As usual, we extract the asymptotic behavior at z 
with the help of the substitution 

VKC) = (-C)-"V(C) • 

This brings Eq. (2) to the following differential equation for /(C): 

v p 



(f"(() + (1 - 2ii/)/'(C) + 



k 2 h 2 



(i-C) 2 (l-C) 



/(C) = o. 



The singularity at C = 1 suggests replacing 

/(C) = (i-C)M0, 

with 



leading to a hypergeometric equation for the function w(Q: 

C(l - ()w"(() + {(1 - 2iv) - [(1 - 2iv) + 25] C} w'(C) - [S 2 ~ 2iv6 + 
The solution is the hypergeometric function 

w(C)=F[S-i( f jL + u),S + i{ f jL-u), l-2iv- C] 



(8) 
(9) 

+oo (C ^ 0) 
(10) 

(11) 

(12) 
(13) 



vp 



W (C) = 0. (14) 



(15) 



up to some normalization factor. For z — > +oo (C — > 0), this function tends to 1, and the 
solution of Eq. (2) contains only the transmitted wave ip(z) ~ exp[i£>3 + W^] satisfying the 
asymptotic condition in Eq. (4) with B ~ 1. For z — > — oo (C — > — oo), we find the asymptotic 
form of the function ip(z) in accordance with the asymptotic condition of Eq. (4) via the 
Kummer transformation of the hypergeometric function w((), with the coefficients 



A 1 = 
A 2 = 



r(l - 2iv) r(-2ifj) 



T[5 - i(fi + v)] T[l - 5 - i{n + v]) ' 
T(l -2iu)T{2ifi) 



(16) 
(17) 



T[5 + - u)]T[l - 5 + - u}) ' 
From these we determine the reflection and transmission coefficients for the Sauter potential: 

1 



R 



l^il 2 ' 



T 



l^il 2 ' 



(18) 
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In our discussion of pair production, we shall focus mostly on the reflection coefficient R. To 
find a simple expression for it, we relate the total kinetic energy Pq^ — Pq + ^ = 2v to the 
parameters v,k of the Sauter potential (3) via Eq. (13). For Av 2 /k 2 h 2 < 1, the square root 

6 = 



! \ — 4v 2 /k 2 h 2 /2 in Eq. (13) is real. After substituting A\ and A2 from Eq. (17) into 
Eq. (18), the reflection coefficient takes the form 

cos 2 n5 + sinh 2 n(fx — v) 



R 



cos 2 7r5 + sinh n(fi + i/) 
where we have used the relation T(l/2 + z)T (1/2 — z) = it / cosnz. 

In the limit k — > 00 where 5 — > 1/2, Eq. (19) reproduces the well-known coefficient [7]: 



(19) 



R 



for the reflection of a relativistic particle off a potential step 

1 la , z > 



<f>(z) =a[6(z) - 9(- 



-z = 



z < 



(20) 



(21) 



k = \^Av 2 /k 2 h 2 - 1 



which arises from Eq. (3) in this limit. 

For Av 2 /k 2 h 2 > 1, the square root 5 in Eq. (13) becomes purely imaginary, 5 = m, where 

(22) 

is real. Then the reflection coefficient in Eq. (18) with A\ and A 2 of Eq. (17) takes the form 

cosh 2 ttk + sinh 2 n(fj, — v) 

R = 7-0 : — r~9 — r , (23) 

cosh ttk + sinh 7r(// + v) 

where we have used the relation r(l/2 + iy)T(\/2 — iy) — n/ coshvry. 

A special case of uniform electric field along the z-direction is included here in the limit 
k — > 0, v — > 00 with t>/c = — eP = const. It corresponds to a linear potential 

4>(z) = -Ez , E < , 

in Eq. (3). In this limit, k ~ v/kh becomes very large, and \x 
so that the reflection coefficient (22) reduces to 



R 



1 + exp 



n(p 2 ± + m 2 ) 
h(e\E\) 



(24) 

v w k/2 - (p^ + m 2 )/4ae|E|, 



(25) 



We have assumed in this section that the particles can pass through repulsive potential 



(-) 



Pn(z) 



barrier with the same (positive) sign of the initial and final energies P c 

Po + v > and Pq + ^ = Po(z)\ z ^ +QO = Po — v > 0, respectively. In the opposite situation, when 

the initial energy Pq ^ = po + v > is positive but the final Pq + ^ = po — v < becomes negative, 
the famous Klein paradox arises. For a particle moving from left to right this means that the 
region of large positive z can only be accessible to antiparticle. Then a non-zero transmission 
coefficient must be present even for a strong potential Av > 2m, where Av = e<p^ —ecj)^ = 2v 
is the potential energy difference at infinity. As has been explained in Refs. [3-8], this happens 
ioi po + v >m and p — v < —m, implying the spontaneous production of particle-antiparticle 
pairs with the total energy difference P^ — Pq^ = 2v > 2m. 
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3 Pair production 



The pair production from the vacuum is now derived as follows. The average number of created 
pairs in the scattering process yields the same result as the imaginary part of the effective action 
in Schwinger calculation [9,10]. This means that the probability for vacuum to remain a vacuum 
under the influence of the external potential, i.e., the vacuum persistence probability, is related 
to the reflection coefficient by 

l(0|0)| 2 = II R P = exp(j2 ^R P ) , (26) 

P0,P± \P0,P± / 

where the product and the sums are taken over all relevant quantum numbers po and pj_ of the 
created particles. Correspondingly, the pair production probability is 



P = 1 - exp ( ^ lntfJ « - lnR p- 

\P0,P_L / P0,P± 



(27) 



In a box-like volume V±T with V± — J d 2 x± = J dx\dx 2 being the area of the potential step 
transverse to the z-direction, and T the total time, the sum over p and p_i_ becomes an integral 

£ h * = W 7^ /?£*) (28) 



P0,P± 



where the reflection coefficient R(po, P±) is defined by Eq. (19) or (23), and p± = |pj_|. The 
rotation invariance around the third axis reduces the integral / d 2 p± to ir J dp\. The remaining 
integral over p and p\ is done over the Klein region 



*o ) = Po + v > \Jp\ + m 2 , P (+) = Po - v < -\fp\ + m 2 . (29) 

The change of the sign of Pq{z) is necessary for a vacuum pair production by the Sauter 
potential (3) with v > m. 

The pair production probability is now completely defined by Eqs. (27)-(29). The probability 
per unit area and unit time is 

1 r(v 2 —m 2 ) rv— y/p^+m 2 

w± = - 3 / dp 2 ± dp lnR(p , p\) , (30) 

2(27t) 2 h JO J-v+y/p 2 ± +m 2 

where the integration region in the (p^_,po)-plane i s shown in Fig. 1. For the actual calculation, 
we interchange the order of integration in Eq. (30) to 



1 

w±= 



2(2tt) 2 ^ 3 



f(po+v) 2 — m 2 rv—m f{po~ v) 2 — m 2 

dPo / dp 2 ± + dp dp\ 

-v+m JO JO JO 



lnR(p ,p 2 ± ). (31) 



We now replace p — > — p in the first integral of Eq. (31) and make use of the symmetry of the 
reflection coefficient R(— Po>P±) — R(Po,P±) under the interchanging /i <-> v in Eqs. (19) and 
(23). The new integration region in the (po,Fj_)-plane is shown in Fig. 2. As a result, we obtain 
for the pair production rate per unit area in the Sauter potential (3) the integral representation 

I rv—m Mpo—v) 2 —m 2 

w± = - 3 / dp / dp\ mi?(p ,pi) ■ (32) 

(27!-) 2 n JO JO 
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v — m 




-v + m 

Figure 1: In the (p 2 L ,po)-pla,ne, the integration covers the positive region restricted by two 
intersecting parabolas Po = v — \Jp\ + rn 2 and po = —v + \Jp\ + rn 2 with horizontal axes of 
symmetry above and below the p^-axis for v > m. 




v — m 



Figure 2: In the (j>o ; pi)-plane, the integration covers the region under the left branch of the 
parabola p\ = (po — v) 2 — m 2 in the first quadrant for v > m. 

4 Potential Step 

For a potential step of the Sauter type (3), the reflection coefficient as a function of po and p\ 
is given by Eq. (20). Its logarithm reads 



\nR(p ,p 2 ± ) = 2 In 



Q+(Po,P±) ~ Q-(Po,P±) 
q+(Po,P±) + q-(po,P±) 



(33) 



where q±(po,p 2 ±) = p^ = J (po ± v) 2 — (p 2 ± + m 2 ) with q + > g_ for p > 0. As in Eq. (8), the 
functions q±(po,p 2 L ) satisfy the constraint 



1 2 r n 2 

q+{po,p±) - q-(po,p±) =4vp . 



The constraint suggests introducing a parameter 9 so that 



q+(Po,P±) = \/4vpo cosh^, q-(p ,p 2 ±) = sf^vpo sinh6». 
This allows us to express Eq. (33) in terms of 9 as 

lni2(po,pl) = -40, 

where 

q+(Po,P±) +q-(po,P±) 



9 = 9 + (p ,p 2 ± ) = \n 



V^upo 



(34) 
(35) 
(36) 

(37) 
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Figure 3: The integration covers the upper region restricted by the logarithmic curve 9 = 
9 + (p ,0) and the positive p -axis in the (p , 0)-plane. Below the p -axis lies the alternative 
region restricted by the mirrored curve 9 = 9_(p Q , 0) = — 9 + (p , 0). 



with < 9 < 9 + (p , 0) for < p\ < (po — v) 2 — m 2 . It is useful to eliminate the variable p\ in 
favor of 9 in the first integral of Eq. (32). 

Alternatively we could have defined 

Q+(Po,P±) - q-(po,pi) 



= MPo, Pi) = In 



V^vpo 



(38) 



with 9_(p ,0) < 9 < for < p\ < (p — v) 2 — m 2 where 9_(p Q ,p 2 ± ) = —9 + (po,p 2 ± ) due to 
Eq. (34). This definition corresponds to Eqs. (35) and (36) with 9 — > —9. 

The change of the measure in the first integral of Eq. (32) due to substituting p\ — ■> 9{jpQ ) p 2 1 ) 

is 



dp\ = —AvpQ sinh2# d9 . 
With Eqs. (36) and (39), the pair production rate (32) takes the form 



w± 



n 2 h 3 JO 



dpoPo 



-(po,o) 



d9A9 sinh2#. 



(39) 



(40) 



where the integration region in the (po, #)-plane is shown in Fig. 3. The integrals in Eq. (40) 
are now straightforward to do. The right-hand integral yields 



Mpo.o) 



d9A9 sinh2# 



1 



2vp c 



(q 2 + +q 2 _)\n 



q+ + q- 



\/4upo 



-Q+ q- 



(41) 



p±=o 



After this, the remaining integral in Eq. (40) becomes a combination of elliptic integrals via 
the substitution po = (v — m)t with < t < 1. The last term in (41) leads to 



rv—m 

/ dp q + (po,0)q^(p ,0) 
Jo 



a+a 2 _ hia) 



where ^{cx) is the dimensionless t-integral 



J 2 (o) = / dtJ(l-t 2 )(l-a 2 t 2 



3a 2 



1 -a 2 )K(a) - (l + a 2 )E(a) 



(42) 



(43) 
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and K(a), E(a) are complete elliptic integrals of the first and second kind, respectively, with 
the argument a = a^/a + < 1, where a± = v ± m. With a little more effort we transform 
Po-integral over the first term in (41) into a ^-integral: 



In 



g+(p ,0) + g_(p ,0) 



= a + a 2 ii(a) , 



(44) 



rv — m ,- 

y o [g+(Po,0) + g!(p ,0) 

where ii(a) is the dimensionless integral 

I^a) = £dt(l + at 2 ) 1 2 In ^(1 + t)(l + at) + ^(1 - t)(l - at) - hit - ln2(l + a)| . (45) 

This can be expressed in terms of the elliptic integrals of the first and second kind F(tp, 1/a) 
and E(ip, 1/a) as follows: 



9a 



1 - a)(4 + a) F(ip, 1/a) + (1 + a(3 + a)) E(ip, 1/a)] , 



(46) 



with ip = arcsin a. 

Finally, collecting all contributions in Eq. (40) yields the pair production rate per area 

2 



v 



w± 



+ 



3tt 2 ^ 3 
a 



m\ 3 
1 + - 

u / 



l-a 2 )K(a)-(l + a 2 )E(a) 



1 - a) (4 + a) F (y>, ^ + (l + 3a + a 2 ) E (y>, J . (47) 



5 Sauter potential 

We employ now Eq. (32) to compute the pair production probability for the Sauter potential (3) 
where the reflection coefficient is defined by Eqs. (19) and (23) for all values of the parameters v 
and k. In order to illustrate the calculation, we specify these, for example, as (2v) 2 /k 2 h 2 > 1. In 
this case, the parameter k defines the inverse width of the electric field, whereas the parameter 
v governs its size \E\ = vk/e, whose maximum is \E C \ = m 2 c 3 /eh ~ 1.3 x 10 18 V/m. The limit 
k — > with fixed vk reproduces the linear potential due to a constant electric field. 

The reflection coefficient of the Sauter potential with Av 2 /k 2 Ti 2 > 1 is given by Eq. (23). An 
equivalent form of this is 

^ cosh ir(/i — v + k) cosh 7r(/x — v — k) 
cosh 7r(/x + v + k) cosh 7r(/i + v — k) ' 

where \i = (J,(po,p 2 ±) and v = v(po,p 2 ±) are the functions of po and p\ defined by Eqs. (5), (6) 

and (7) with the constraint (8), while k = \J Av 2 — k 2 h 2 /2kh is a constant. Taking logarithms 
of Eq. (48) leads to the expansion 

ln R = — 4 ^2 cos h 2rrnK sinh 2nnfi sinh 2imv . (49) 

n=l n 

The right hand side is found by replacing each logarithm of the hyperbolic functions as 
ln(2coshx) = x + ln(l + e~ 2x ) = x + Y^=i(— l) n+1 e~~ 2nx / n, and combining all sums. 
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With Eq. (49), the pair production rate per area (32) takes the form 

1 00 (—i) n + 1 

w ± = y ± '- cosh(2vrn/t) J (n) , (50) 

where are the integrals 

rv-m rv~m r(p -v) 2 -m 2 

J (n> = / dpo I {Po) = / dp / sinh27rn/i(po,Pj_) sinh27rnz/(po,Pj_) • (51) 

./() ./n JO 

The region of integration is shown in Fig. 2. 

A physically more instructive quantity than the production rate (32) can be obtained by 
dividing w± by the width of the potential step. For the Sauter potential the width is defined 
by the ratio 

/oo 
dzE 2 (z)/E 2 max = 4/3k. (52) 

Thus we obtain the pair creation rate per volume of nonzero field 

w = w±/L = (3k/4)w±. (53) 

We perform the ^-integration in Eq. (51) by changing, for each n separately, from the 
variable p\ to the dimensionless one 9 defined by 

p2 _>0= [2vmz/(p ,pi)] 2 (54) 

with < 9 < {n) {p o ), where 

W (po) = [2nnu(p , 0)] 2 = (nn/kh) 2 [(p - vf - m 2 ] . (55) 

For each n, Eq. (54) yields 2irnis(p ,p 2 L ) = \J~9~ and dp\ = — (kh/nir) 2 d9. Let us also introduce 
the dimensionless functions 

9^(p ) = (2nn/kh) 2 vp , (56) 



so that 27infi(po,p 2 L ) = yO^ipo) + 6 due to Eq. (8). Then ^-integrals take the form 

f(P0—v) 2 —m 2 

I (Po) = / dp 2 ± sinh 27rn/i(po, P±) sinh27rnz/(po,pl) 

Jo 

= 7^^J° (P0) ^{cosh^(«)(p ) + 9 + y/6j- cosh^M^) + 9 - y/fty. (57) 
The two terms in Eq. (57) can now be combined into a single integral as follows. We 



substitute t = \/9( n \p ) + 9 + \J~9~ with ^ (n) (Po) < t < 9+\po) in the first term, and 
t = ^9^(p ) + 9 - V9 with 9 { l l \p ) < t < ^9( n \p ) in the second, where 

9t\p Q ) = (27^)^(^,0) ± !/(po,0)]. (58) 

Noting that \/# = ±\t/2 — 9^ n \p )/2t] in the first and the second substitution, respectively, 
while d9 = (t/2){l - [9 in \p )] 2 /t 4 } dt in both cases, we obtain 
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Evaluating this integral yields 



/ (n) (Po)=/| n) (Po)+/i n) (Po), 



(n), 



with 



/ kh\ 



± 



4 n) sinhC J -cosh^' 



(n) 



(n) e (n) 



Chi#i n) - 



sinh4 n) cosh^i n) 



9 (n) 



(60) 



(61) 



where 9±^ is short for 9±\p ), Chi6^ are the hyperbolic cosine integrals, and the last two 
terms represent the leading terms in their asymptotic expansions for large arguments 9±\ 

Having obtained I {n) (p ) = I^(0 { + \p ) , 9 { " \p )), we are left in Eq. (51) with the sum 



jw= r m dp i^( P0 )=4 i) +j^, 

Jo 

of the integrals over the rather lengthy functions (61): 

, . rv-m , . 

4 n) = / dp,it\p ). 

Jo 



(62) 



(63) 



However, this sum can be combined into a single integral by subjecting J± ^ in Eq. (63) to a 
change of variables po — > £ provided that we define the new dimensionless integration variable 
£ as follows. 



In order to transform the integral J+ \ we introduce the dimensionless variable 

C(Po) = MPo,0) + i/(p ,0) , 



(64) 



with yv(v — m) /%k < £ < \Jv 2 — m 2 /hk for < po < — m ) an d v > m, where /i(po? 0) = 

\J (Po + ^) 2 — m 2 /2hk and ^(po,0) = \J (p — v) 2 — m 2 /2hk. Noting that m 2 < [v 2 — {hkQ 2 } 
within these limits, we resolve Eq. (64) in terms of p as 

1 1/2 



m 



[v 2 - (hkO 2 } 



(65) 



with positive £ due to po > 0. To determine 9±\po) of Eq. (58) in terms of £ by means of 
Eq. (65), we find first for a given £ the positive square roots 



\J (po ± v ) 2 — m 2 = hkl; ± v 1 1 
Combining these yields 



m 



1/2 



b 2 - (hkO 2 } 



> o. 



//(p , 0) + is(p , 0) = £ , /i(p , 0) - v(po, 0) = — I 1 



1/2 



[w 2 - (^o s 



(66) 



(67) 



From Eq. (67) we obtain, finally, the functions 0±\£) = #±^(po(0) to be substituted instead 
of 9±^ (p ) into the first integral jl™- 1 as 

s 1/2 



^(0 = 2^, ^(0 



27rn 



■t; 1- 



[v 2 - (/ifcf) 5 



(68) 
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Note that the inequality m 2 < [v 2 — (Tik^) 2 } ensures the positivity of expressions under the 
square roots in Eqs. (65)-(68). 

In order to treat the second integral in Eq. (63), we define a new integration variable £ 
similar to Eq. (64): 

£(Po)=MPo,0)-z/(po,0), (69) 

with < £ < \Jv{y — m)/hk for < po < (y — m) and v > m, where /x(p ,0) = 

\J (Po + v ) 2 — m 2 /2Tik and v(p o ,0) = ^/ (Po — v) 2 — m 2 /2hk. Since m 2 < vm < [v 2 — (fr&£) 2 ] 
in these limits, we solve Eq. (69) in terms of p in the same way as in Eq. (65), leading again 
to 



m 2 



1/2 



M0 = m[i- [v2 _ {hkm i , (to) 

with positive £ due to po > 0. To determine 0±^(po) i n terms of £ by means of Eq. (70), we 
find now for a given £ the positive square roots 



(71) 



It follows from Eq. (71) that 



1/2 



MPo,0)-KPo,0) = £, //(Po, 0) + u(po, 0) - ^ |l - [y2 _ m m)2] } (72) 
This yields ^i - ^) which replaces 0±\po) in the second integral as follows: 

^fH^F^blf • ^««. P3) 

Again, the inequality m 2 < [t> 2 — {Tik^) 2 } ensures the positiveness of the expressions under the 
square roots in Eqs. (70)-(73). 

We now go from the integration variable po to £ in Eq. (63) for and separately. 
Substituting these in Eq. (62) yields 

v 2 1 dpo(0rf(0 



I d£ ^ ^ [cosh (27m + 2vr< sinh (2™$)] 

ft 2 fc 2 « 
(27rn) 2 io * rf£ 



+ 7^ %P t cosh ( 2 ™ " 27r < sinh ( 2 ™ 01 , (74) 
[limY Jo at 



where 

f = Vf 2 - m 2 /^, (75) 
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and the function po(£) is given by Eq. (65) [or (70)]. Equation (74) can be simplified by partial 
integration thanks to the vanishing of the function po(0 a t the endpoints of integration. This 
yields 



jW = (hk) 3 cosh (27m 0, 

Jo 

where the function /(£) has the form 

/(0=ePo(0/^-^ 2 Po(0/3W 5 e 3 . 
After inserting po(0 from Eq. (65) [or (70)], this reads explicitly 



/(0 = 



e - e 



e-t 2 + (m/hky 



+ (m/hky 



e-e 



e - e + {m/hky 



3/2 



(76) 



(77) 



/(-*)■ (78) 



Note that the integral over this function vanishes: fjjd^ /(£) = 0, so that = 0. 
Alternatively, we may introduce the variable rj = £/£ to arrive at the integral 



_ ^j v 2 _ m 2 y dfjgfjj) cosh (2irn£ r^J , 



with 



= /(to))/? = ^ 



1 - 77 2 + 7 2 



1/2 



1+7 2 



-i 3/2 



1 - T] 2 + 7 2 



where 7 = m/hk£ = mj ' \Jv 2 — m 2 . The integrals are all functions of v,m, and fc. 



(79) 



(80) 



6 Pair Production Rates 

Let us first check our final expression (50) with (79) by going to the limit of a constant electric 
field k — > 0, v — > 00 with vk = e\E\, where the exact result is known. In this limit, the 
parameter 7 becomes small and can be neglected, and the integrals (79) become approximately 



Jin) 



sinh(27r£n) cosh(27r£n) sinh(27r£n) 



(27r£)n 



(27r0 2 n 2 



(81) 



Inserting these into Eq. (50) leads to the following hyperbolic sums 



E 

n=l 



"I) 



n+1 



cosh(nA±) 



-^Li„(-e A ±) =p ^(-e- A ±) , 



(82) 



sinh(nA±) 

with v — 2, 3, 4, where A± = 27r(« ± £) and LL,(z) are the polylogarithm functions 

00 n 

Li ^) = E— . ^ = 2,3,4. 



71=1 



(83) 
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Note that the constant field limit corresponds to large arguments in Eq. (82), since k ~ v/hk, 
£ ~v/hk — m 2 /2vhk ~ k — p/27r, where 



p = nm 2 /he\E\ 



(84) 



We must exploit therefore the analytic continuation of the polylogarithm functions defined by 
the series (83) into the region \z\ > 1. By taking advantage of the formula [15]: 



LU-z)+e™U v ( 



I» H 2 2ttz 



where ((is, q) is the Hurwitz zeta function 



(85) 



CM) = E 



n= (* + ?)' 

we bring the right-hand side of (82) to the form 



(86) 



1 (2*)" 



C 1 - V, 



1 A. 



2I» 

For v — 2, 3, 4 this reads explicitly, 

oo / j\n+l 

X) — ~ — sinh(nA ±y 

n=l 



2ni 



^ +e iw Li I/ (-e- x± ). 



(87) 



E 

n=l 

oo 

E 

n=l 





n 2 


(- 








(- 




n 4 



cosh(nA±) 
sinh(nA±) 



7T 

12 

^2 



+ ^ + Li 2 (-e- A ±), 



^A ± + ^±-Li 3 (-e- A ±), 
12 ± 12 ^ ; ' 



(88) 



77T 4 7T 



A' 1 



\ 2 + — 

720 ' 24' ± 48 



+ — A± + + Li 4 (-e 



)• 



Substituting Eq. (50) with Eqs. (81) and (88) into Eq. (53), we obtain the approximate pair 
production rate per nonzero field volume 



w = w±/L 



1 



3A; v 3 \ 2tt 2 . - 1 . _> . 1 . \ . 



+ 



(2*0 
1 

(2<^ 



9 -3i-e- A+ )--Li 3 (-e- 



7T 27T / 9 - 9 \ 1 

_« + __ /c ( /c 2 + 3^__ Ll; 
-(<) + — (<)( K 2 + £ 2 ) + -Li 4 (-e 



Li 4 (— e~ 

2 v 



(89) 



We now take the constant-field limit k — > 0, v — > oo at t> /c = e|i?| fixed, where A_ = p remains 
finite, while A + = 47r/t — p ~ 47r/t — > Arrv/hk ~ 4ireE/k 2 h tends to infinity, so that Li J/ (— e~ A+ ) 
with z/ = 2,3,4 vanishes. Moreover, the polylogarithm functions Li„(— e~ A ~) with z/ = 3,4 do 
not contribute because of vanishing prefactors. All divergent terms cancel each other. Thus we 
obtain the pair production rate per nonzero field volume, which for constant field is the total 
volume: 



w 



w ± 3 (e|£|) 2 

~L ~* ~4 6n 3 h 2 c 



U 2 (-e-") 



(90) 
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Here we have inserted L = 4/3/c from Eq. (52). The division by L is essential for getting a finite 
result in the constant-field limit. For completeness, we have reinserted in the final expression 
the light velocity c to verify the complete agreement with the result of Heisenberg and Euler [1], 
Schwinger [4], Nikishov [10], and many others (see e.g. [13] and references therein). 

For arbitrary k, the integral (76) cannot be evaluated in closed analytic form. In order to 
obtain an approximate rate formula we insert from Eq. (76) into Eq. (50) and interchange 
the order of summation and integration to find the expansion 

L.3 r ( oo (_-i\n+l 

™.L=o-2 / £ {cosh [27m(£-«)]+ cosh [27m(£ + «)]} . (91) 

Jo n=l U 

The integral is simplified with the help of the summation formula 

V - - cosh(n:r) = - fin f 1 + e - *) + In (1 + e x )] , (92) 

„=i n ' J 

which permits us to bring the general rate to the form 

W± = 2^L^ { M I 1 + e ~ Mi ~ K) ] + ln I 1 + e M6+K) ] } , (93) 

or, by the symmetry of the function f{£,) = /(—£), to the more symmetric form 

k 3 r€ r i 

W±= 2^ J-^ m M I 1 + e ^~ K) ] ■ ( 94 ) 

Finally, integrating this by parts, we find 

k 3 ft 1 



9(0 = ~ ~ ' I —2 = ~9(-0 ■ (96) 



where the function g(£) vanishes on both ends. Explicitly, it reads 

i (?-^ 3/2 

(e - e + m 2 /h 2 k 2 ) 

Remarkably, the second function under the integral (95) resembles a Fermi distribution. Indeed, 
we are going to show that the calculation of the integral (95) can be done by a method familiar 
to low-temperature expansions in statistical physics [16]. The condition necessary for pair 
production v = e\E\/k > m implies that the parameter k lies in the interval < k < em/%, 
where 

e=\E\/\E c \ (97) 

and E c = m 2 /eh (in natural units with c = 1) is the so-called critical field for which the work 
over two Compton wavelengths 2H/m can produce the energy 2m of a pair. At the upper end 
k = em/h of the above interval, the rate (95) vanishes, since £ becomes zero. For the calculation 
of the exact pair production rate (53) from Eq. (94), we introduce the dimensionless parameter 
k = hk/em, where < k < 1, to rewrite Eq. (95) in terms of the dimensionless variable 
£ = £hk/v. This brings the production rate (53) to the form 

W = " (iji I-i di ~ 9{i) eM*-*)U + 1 ' (98) 
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where 

£ = yf\ - ~k 2 , k = ^l-{e~k 2 /2) 2 , (99) 
and the dimensionless function g(£) reads 

f (.?- - ; 2 ) 3/2 
(i-e) 

We expand this function into a power series 



~9(i) = I , z A/2 • ( 10 °) 



sf(l) = Ew(6f n+1 » ( 101 ) 

n=0 

with the coefficients 

fli(S = -|(3-?)/2, »5(e) = 3(l-?) 2 /8|, •••• (102) 

Substituting the expansion (101) back into Eq. (98), we encounter the odd-moment integrals 
of the Fermi distribution 



?2n+l 

These can all be found exactly. Performing the integrals yields a binomial expansion 



/ ^Wi ' n -°- (103) 



M 2n+1 (0-^ (-1) (2vr )"T(2n-m + 3) Llm(6 " (1 ° 4) 

Here Li m (£) are the linear combinations of the polylogarithm functions 

n m (|) = Li m (-AW(0) + (-I)™-* U m (-X-\i)) , (105) 

with arguments 

A«(0 = e^), \<-\i) = X + \-i)=S-t>, (106) 

where 

p(i) = 27r(i-R)/e~k 2 . (107) 
The exact production rate of Eq. (98) becomes now the sum 

(e\E\) 2 1 



w 



J292n+i(OM 2n+1 (0- (108) 



Anh 2 (eP )2 ^ 

By making use of Eqs. (104), we rewrite this as an expansion over the polylogarithm functions 

w = - YJ- £ Ul) (^T" 2 n m (0 , (109) 
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where the coefficients c m (|) are polynomials of £ 

^«»-| (2^rpn (2 :^3) ^-«»^" 1+2 - < 110 ' 
Together with Eq. (101) these read explicitly, 

Cl (|) = ( 5 e 4 /2567r)(-l-4f/5-6f/5-4e 6 + 7e 8 + •••), 

C2(|) = d 3 /5127r 2 )(125 + 84^ 2 + 102^ 4 + 260^ 6 -315f + •••), 

c 3 (|) = (3£ 2 /1287r 3 )(15- 10| 2 — 32| 4 — 110| 6 + 105| 8 + ---), (111) 

c 4 (|) = (3^/2567r 4 )(-205 - 70f + 132£ 4 + 910£ 6 - 735£ 8 + • • •) , 

c 5 (i) = (1/256tt 5 )(1347 + 1296£ 2 - 18£ 4 - 9240£ 6 + 6615£ 8 + •••), 



The series expansion given by Eq. (109) converges well for small k. Here the parameter 
/?(£) in Eq. (107) becomes ~ — n/e + n(ek 2 )/A + ■ ■ •, where the first term is equal to — p 
of Eq. (84), and the polylogarithm functions Li m (— A^(£)) with m > 2 will be suppressed 
by powers of k. The parameter p~(—£) = — 27r(£ + k)/ek 2 tends to minus infinity, so that the 
polylogarithm functions Li m (— /V ■*(£)) for all m yield the exponentially small contributions. 
By means of Eq. (99), the coefficients (111) are the polynomials of small k 

Cl (k) = {l/n){-~k 2 /4 + 5~k 4 /4-43~k e /16 + ---), 

c 2 (~k) = (1/tt 2 )(1/2-3P/8-27A;78 + 589fc 6 /64 + •••), 

c 3 (k) = (1/tt 3 )(-3/4 + 3P/8 + 189A;716- 1011^/32 + •••), (112) 

c 4 (k) = (l/7r 4 )(3/8 - 2655fc 4 /64 + 6555^/64 + •••), 

c 5 (jfe) = (l/7r 5 )(279fc 4 /2- 19155^/64 + •••), 



Finally, this yields the probability rate (109) as a series expansion in powers of small k 



w = 



(e\E\f 



471-r 



' Li 1 (-e^) + - i T Li 2 (-e-0 



47TC 



2tt 2 



-Lii(-e-0-^Li 2 (-e^) 



Aire 



3e 

47T 3 



Lia(-e^) 



+ 



■}• 



(113) 



where the leading term is already an excellent approximation. Note the coincidence of the 
second term in the first brackets with the probability rate (90) for a constant-field limit k — > 0. 



7 Conclusion 

We have calculated an exact expression for the production rate of charged scalar particle- 
antiparticle pairs from the vacuum by the Sauter potential. For an arbitrary potential barrier, 
the rate was related to the scattering amplitude on the barrier, and expressed as an energy- 
momentum integral over the logarithm of the reflection coefficient. For the Sauter potential, we 
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have evaluated this integral and checked the result by recovering the known limits of a sharp 
step potential and of a uniform electric field. 
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